SUBSTITUTE SPECIFICATION 



"Method for determining reflections in an area". 
* * * 

DESCRIPTION 

BACKGROUND OF THE INVENTION 

The present invention refers to a method for determining reflections in 
complex environments, by means of ray tracing. 

The rapid tracing of rays, that follow reflection laws similar to geometric 
optics laws, is a problem of considerable interest in numerous fields of 
application. In particular, the techniques capable of carrying out this operation 
find immediate application in the simulation of acoustic reverberations in 
complex environments, in the photo-realistic rendering of models of three- 
dimensional scenes, or in the prediction of the electromagnetic propagation of 
wireless systems in the presence of multiple reflections. 

The acoustic or electromagnetic path curves follow paths that, with 
excellent approximation, meet the geometric optics laws. According to these 
laws, a ray that has an incidence on a planar reflector with angle a in relation to 
the normal to the plane is reflected by this in coplanar direction with the plane 
formed by the incidence vector and the vector normal to the reflector, which 
results specular in relation to the latter vector (output angle equal to a -a). A ray 
can undergo many reflections by the walls present in the environment, each 
time diminishing until it becomes of negligible entity. Given a source and a 
receiver however located in an environment with various reflectors, many paths 
can exist (in direct or indirect visibility) that connect these points. 

There are various approaches to the problem of tracing rays, such as the 
image sources method, the ray tracing method, and the beam tracing method. 

The image sources method calculates the paths corresponding to specular 
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reflections considering the virtual sources obtained through the specular 
reflection of the sources in relation to all the polygons that represent the 
reflectors of the environment. For each virtual source a specular reflection path 
can be constructed through the iterative intersection of a rectilinear segment that 
links source with receiver through the existing reflection planes. The specular 
reflection paths can be calculated up to an arbitrary order through the recursive 
generation of virtual sources. The main advantage of the image method is its 
solidity, as it guarantees that all the required paths are determined up to a pre- 
assigned reflection order. Against this, the cost of calculating the method grows 
in exponential proportion of the reflection order and the control of the visibility 
conditions (to be carried out for each of the virtual sources) results to be of 
limited complexity only for environments that are not too complex (in fact 
conditions of occlusion, plausible reflection, etc. have to be considered). 

The ray tracing methods are based on the direct tracing of the ray starting 
from the source until it reaches the receiver or vice versa. Whatever point the 
ray tracing starts from, these methods typically require the exploration of a 
finite, but high, number of output directions, to verify that a path exists that 
complies with the geometric optics laws and that reaches the arrival point after a 
certain number of reflections. However, this approach often results to be 
burdensome to calculate as in order to be certain not to excessively undersample 
the space of the possible directions to explore (and erroneously miss the arrival 
point), the number of tracing attempts has to be very high. 

Beam tracing is capable of preventing this problem through the tracing of 
compact beams of rays that originated from the same source (real or virtual). 
This eliminates the risk of missing the arrival point in the tracing operation. 

When a beam has incidence on a reflector, all the rays that compose it are 
reflected according to the geometric optics laws forming a reflected beam that 



results in coinciding with the specular version (in relation to the surface of the 
reflector) of the beam that would be obtained if the reflector was not present. 

In general, however, the beam has incidence only in part on the reflector, 
with the result that only a portion thereof is reflected while the remaining part 
proceeds undisturbed in its original path. In this case the beam separates into 
two or more sub-rays in correspondence with the reflector. To carry out beam 
tracing thus means to construct the tree of the reflections and of their branches 
that the beam undergoes in its path in space. Clearly, the ray undergoes a natural 
attenuation that grows with the increase of the length of the path followed. A 
further attenuation is introduced by the reflection walls. This enables the growth 
of the tree to be interrupted where the attenuation of the beams reaches a preset 
threshold. 

The reflection tree or beam tree is constructed so as to gather in gerarchic 
mode all the information needed for the reconstruction of the acoustic paths 
existing between the source and an arbitrary point. In particular, for each branch 
point of the tree the position of the source (real or virtual) that generates the 
reflection to that branch point has to be specified, the geometrical parameters 
that characterise the extension of the beam, and the information need to identify 
the reflector - if any - on which the beam will have incidence in following its 
path. It is assumed, without loss of generality, that the beam has been 
preliminarily divided into adjoining sub-beams. 

Once the beam tree associated to a given position of the source has been 
constructed, the determination of all the beams that hit a certain point of the 
space (receiver) is immediate. Going back up the tree it is then possible to 
derive the direct or reflected paths that link it to the source. This operation 
results, in fact, in being very rapid as it only requires consultation of the data 
structure used to memorise the beam tree. 



It can be seen however that a movement of the source or even of only one 
of the reflectors entails the recalculation of the tree, an operation that is 
generally more burdensome than its simple consultation. 

SUMMARY OF THE INVENTION 

In view of the state of the technique described, an object of the present 
invention is to speed up the construction of the beam tree when the source 
moves, through the separation of the phases of analysis of the visibility and 
tracing of the beams. 

In its first aspect the present invention refers to a method for determining 
the region of visibility between at least a first reflector and a second reflector 
comprising the following phases: representing said first and second reflector in 
a system of coordinates (x, y, z); said method is characterised in that it further 
comprises the phases of: carrying out a similar transformation of said system of 
coordinates (x, y, z); determining the region of visibility of said second reflector 
in relation to said first reflector as the set of the parameters of the straight lines 
that link a generic point of said first reflector with a generic point of said 
second reflector. 

Preferably, the affine transformation comprises the execution of a rotatory 
translation of said system of coordinates (x, y, z) so that said first reflector is 
placed on the plane x = 0. 

Preferably, the affine transformation comprises the execution of a scaling 
down of said system of coordinates (x, y, z) so that said first reflector assumes a 
preset dimension. 

Advantageously, the method comprises a third reflector; and the following 
phases: determining the region of visibility of said second and third reflector 
seen from said first reflector; in the case that there are overlapping regions of 
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visibility: tracing a semistraight line in said system of coordinates that has 
parameters comprised in the overlapping region starting from said first reflector; 
determining what is the succession of the intersections between said second and 
third reflector; assigning the overlapping portion of the region of visibility to the 
5 reflector having priority succession order; said regions of visibility represent a 

visibility diagram. 

In particular the method is characterised in that it represents said first 
reflector and second reflector in a system of coordinates (x 5 y) by means of 
segments; carrying out an affine transformation of said system of coordinates 

10 (x, y) capable of bringing said first reflector to assuming coordinates of the 

extremes in the points (0, m) and (0, n) with m < n; representing said second 
reflector by means of the following system of equations x = et + f 5 y = gt + h 5 
with 0 < t < 1 ; representing a generic straight line by means of the parameters a, 
b of the equation y = a x + b; determining the region of visibility by determining 

1 5 all the straight lines that pass through a generic point of the first and of the 

second reflector combining the previous equations and obtaining the following 
system gt + h = a(et + f) + b, with 0 < t < 1 and with m < b < n. 

In its second aspect the present invention refers to a method for determining 
a beam tree of the beams of the rays on a plurality of reflectors that comprises 

20 the phases of prior determination of the set of the regions of visibility from each 

reflector in accordance with the previously described method. 

Advantageously, the method for determining a beam tree of the beams of 
the rays on a plurality of reflectors comprises the phases of prior determination 
of the visibility diagram in accordance with the previously described method; 

25 positioning a source in a system of coordinates (x, y); determining the areas of 

the reflectors illuminated by said source in said system of coordinates (x, y); 
memorising the coordinates of said areas of the reflectors illuminated by said 
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source; representing a virtual source for each reflector illuminated; applying 
said affine transformation of said system of coordinates (x, y) to said virtual 
source and to said illuminated region; said transformed illuminated region is 
represented by means of a segment of extremes (0, bo) and (0, bi); representing 
said transformed virtual source in the space of the parameters (a, b) by means of 
the equation y = a x + b, and the illuminated region by means of the disequation 
bo < b < bi; the system of equations y = a x + b and b 0 < b < bi represents a 
segment of illumination in the space of the parameters (a, b); intersecting said 
segment of illumination with said visibility diagram obtaining sub-segments and 
thus sub-intervals of the interval bo < b < bi; said sub-intervals will represent 
portions of the illuminated region and each one will illuminate a new reflector. 

In its third aspect the present invention refers to a program for a computer 
comprising a program code that carries out all the phases of the method for 
determining the region of visibility between at least a first reflector and/or of the 
method for determining a beam tree when the program is carried out on a 
computer. 

In its fourth aspect the present invention refers to a program for computer 
recorded on a support that can be used by a computer for controlling the 
execution of all the phases of the method for determining the region of visibility 
between at least a first reflector and/or of the method for determining a beam 
tree. 

The present invention permits to carry out the bulkiest calculation 
operations of assessment of the visibility in a preliminary phase that considers 
only the reciprocal positions of the reflectors, independently from the position 
of the source. The visibility is in fact redefined in a suitable parametric space, 
where the information for its assessment can be organised under the form of 
diagrams. The construction of the beam tree is then carried out each time the 



source moves, using the precalculated information of visibility. 

BRIEF DESCRIPTION OF THE DRAWINGS 

The characteristics and advantages of the present invention will appear 
evident from the following detailed description of an embodiment thereof, 
illustrated as a non-limiting example in the enclosed drawings, in which: 

Figures la, lb and lc show an example of rotatory-translation, of 
normalization and transformation from a space with coordinates (x, y) to a 
space with parameters (a, b); 

Figures 2a and 2b show an example of construction of the visibility 
diagram; 

Figures 3a and 3b show an example of construction of a level of the beam 

tree. 

DETAILED DESCRIPTION OF THE INVENTION 

As the method herein proposed refers to the tracing of paths that follow the 
geometric optics laws, phenomenon such as diffraction will not be considered, 
which occurs when the dimensions of the reflectors become comparable with 
the wavelength of the wave that propagates in the environment. This 
assumption, even though correct in the case of high frequency electromagnetic 
waves (light or microwaves), results only approximately satisfied in the case of 
the propagation of acoustic waves, and in fact it has been calculated that these 
models are valid for simulations above 150Hz (corresponding to the wavelength 
of approximately 2 metres, a length that starts to be comparable with the 
dimensions of the walls of a living or working environment). In addition, even 
though it is true that sufficiently wide walls reflect in an almost specular 
manner, rough edges of rigid reflectors tend to generate a diffraction that is not 



negligible. In particular, small objects (a vase, a cup, etc.) tend to have 
significant effects on the acoustic field for sounds whose frequency is above 
4kHz. These objects, however, are not generally considered in the modulation 
of acoustic environments. 

For simplicity of description, hereinafter the phenomenon of the 
transmission through the walls will also be neglected. Hereinafter it will be 
illustrated how it is possible to take account of it through simple modifications 
of the method illustrated herein. 

The entities that are now considered are the sources, the receiver and the 
reflectors. 

It is assumed that the sources and receivers are pointed elements, while the 
reflection walls are described by segments in the two-dimensional case and by 
portions of polygonal planes, for example triangles, in the three-dimensional 
case. Each wall of reflection defines two reflectors, one for each side. This 
distinction, as will be seen, is substantial in that it has consequences on the 
choice of parameters that characterise the visibility. 

It can be said that a source "illuminates" a point if there is a direct path that 
links them. The term "virtual source" refers to the specular reflection of a real 
source in relation to a reflecting wall. In addition, that region of a reflector that 
is directly illuminated by the source considered is defined as "active". The 
active portion of a reflector can be composed of one or more active segments 
(also disconnected because of possible partial occlusions). Each of the active 
segments generates a new reflected beam. 

Identifying the reflectors present in the environment with a whole index i, 
the visibility is defined as a function that associates the index of the visible 
reflector to a view point and to a view direction. This function is thus defined in 
the space of the parameters that characterise view point and direction, and is 
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constant at times. The visibility from a reflector is thus defined as the visibility 
function from all points of the reflector itself. 

The assessment of the visibility of an environment, that is of a set of 
reflectors, consists in the construction of the visibility functions associated to all 
the reflectors of the environment. As it is assumed that the reflecting walls 
reflect on both sides, the visibility of an environment composed of M walls will 
be characterised by 2M visibility functions. 

The preliminary assessment procedure of the visibility of the environment 
through the construction of the visibility by the generic reflector is now 
illustrated. To be able to conveniently define the space of the parameters on 
which the visibility function assumes values, it is preferable to carry out an 
affine transformation of the geometric space that leads the reflector to assume 
preselected position, orientation and, finally, shape. 

As previously illustrated, the visibility function assumes values on the 
space of the parameters that characterise a view point and direction. If the point 
lies on the reflector, the space of the parameters will be four-dimensional in the 
case of three-dimensional environments, or two-dimensional in the case of two- 
dimensional environments. It is assumed that the geometric space has been 
preliminarily submitted to affine transformation that brings the reflector to lie 
on plane x = 0, with reflecting face turned towards the semi-axis x > 0. A 
convenient parametizing for view point and direction is thus given by the 
system of equations 

y =ax+b 

z=cx+d 

where the couple of parameters (b, d) identifies the view point on the plane 
(y, z) of the reflector, while the couple of parameters (a, c) defines the view 
direction. In the case of two-dimensional environments only the first of the two 



equations has to be considered, and the parameters a and b will identify the 
view direction and point respectively. 

For reasons of simplicity of description, only the case of two-dimensional 
environments is now considered. In this case, the affine transformation of the 
geometric space can be defined simply as a rotatory-translation with scaling 
down that transforms the reflector in the segment that links the points(0, -1) and 
(0, 1), with reflection face turned towards the semi-axis x >0. This enables the 
space of the parameters at the interval -1 < b < 1 to be limited. The visibility 
function will thus be a constant function at times in the space (a, b), where the 
values identify the visible reflectors. 

Two reflectors r 0 and ri are considered, having as references respectively 
qo, po and qi, pi, located casually in a space with coordinates (x, y), as shown in 
Figure la. 

The reflectors are roto-translated so as to place the reflector ro in question 
to coincide with the vertical axis y (x=0). Again for practicality, the reflectors 
are normalized, that is the segment having as references p'o and q'o of 
coordinates (0, 1) and (0,-1) has been chosen as representation of the reflector 
ro. Thus a representation of the reflectors as shown in Figure lb is obtained. 
Representing the reflectors in the space of the parameters (a, b) the diagram in 
Figure 1c is obtained. 

In the space of the parameters (a, b), in Figure lc, the space o comprised 
between the four straight lines p"o, p"i, q'o, q"i represents the visibility of r'i 
from r f o. 

The region of visibility of a generic reflector in relation to the reflector of 
reference is defined as that region of the space of the parameters (a, b) to which 
view points and directions from which that reflector is visible correspond. As 
occlusions by other reflectors can occur, this region can result as being empty, 
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or composed of one or more convex polygons. In each region of visibility the 
visibility function will obviously result to be constant. Vice versa, the region of 
visibility associated to the umpteenth reflector, corresponds to the region in 
which the function assumes value i. It is in addition useful to define the region 
of potential visibility of a reflector in comparison to the reference reflector, like 
the region of visibility calculated in absence of all the other reflectors. The 
effective region of visibility will be justly contained in the region of potential 
visibility, as the regions of potential visibility of all the reflectors can overlap. 

The construction of the visibility function associated to a reference reflector 
is carried out in two successive phases: first the regions of potential visibility of 
the single reflectors are assessed, then the regions of effective visibility are 
assessed through determinations of the reciprocal occlusions. 

For determining the region of potential visibility of the generic reflector 
compared to the reference reflector, the parameters of all the straight lines that 
start from any point of the reference reflector (in positive direction) and have 
incidence on the reflecting side of the reflector considered at any point are 
determined. To this object, first of all the generic reflector is defined with the 
couple of equations: 

x=et+f 

y=gt+h 

where t is a parameter comprised between 0 and 1 . Replacing these 
expressions in the equation y = ax + b a region in the space of the parameters is 
obtained, which is characterised by: 

gt + h = a(et + f) + b, 0 < t < 1 

-1 <b< 1 

This region results to be thus the intersection between a beam of semi- 
straight lines and the reference strip -1 < b < 1. 
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When the areas of potential visibility of two reflectors overlap each other in 
the space of the parameters, it means that a reflector partially blocks the other in 
relation to the generic view point on the reference reflector. The second phase 
of the procedure for assessing the visibility function consists in assigning the 
correct index to the regions in which the regions of potential visibility overlap 
each other. This operation can be carried out in various ways. The most simple 
method consists of determining which reflector occludes all the others by means 
of tracing a sample view ray chosen inside the overlapping region. At the end of 
this operation a partitioning of the space of the parameters in areas of effective 
visibility is obtained, to each of which is associated at the most a visible 
reflector. This geometric representation of the visibility function, is referred to 
as a visibility diagram. 

Now four reflectors rj-r4 are, and the visibility diagram of the reflector is 
to be determined. Thus, the reflector 14 is roto-translated and normalized, and a 
representation in space with coordinate (x, y) is obtained as shown in Figure 2a, 
and in the space of the parameters (a, b) as shown in Figure 2b. A source S 
represented in the space of the parameters (a, b) with the straight line S f is also 
considered. It can be seen that in Figure la, the reflector r3 partially occludes the 
reflector rj. An occlusion in the space of the parameters (a, b) is represented by 
an overlapping of areas. Analysing the slopes of the straight lines that enclose 
the overlapped areas it is possible to determine in the correct way the relative 
visibility. In particular, the greater slope dominates the minor slope. 

This operation of construction of the functions (diagrams) of visibility is 
carried out on all the reflectors present in the scene (two per wall) in a 
preliminary analysis phase. 

At this point the real tracing of the beams for the construction of the beam 
tree can be started, an operation that is possible when the position of the source 
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is known. As already partially anticipated, each branch point of the beam tree is 
associated to a beam, and is characterised by the following information: position 
of the real or virtual source, index and illuminated portion of the reflector on 
which the beam has incidence. 

The present method enables in particular the construction of the upper 
branching of the beam tree at the first level to be sped up. The construction of 
the branching of the root of the tree (primary beams) can come about, as in the 
classic method, tracing a beam at 360 degrees from the source and verifying the 
occlusions in the geometric space. 

The beam tracing procedure comes about in recursive mode. After having 
divided the beam emitted by the real source into primary beams, each associated 
to a reflector in visibility, they are reflected and the recursive determination of 
the successive branching is made. 

Given a beam reflected by the umpteenth reflector, the information of the 
visibility diagram associated to it is used to determine the subdivision of the 
beam into sub-beams, one for each reflector in visibility. 

In practice, the procedure is as follows: first of all the transformation used 
to bring the reflector associated to the beam in canonic reference position is 
applied to the virtual source, that brings it to lie on the left of the reflector. At 
this point the interval bo < b < bi that characterises the portion of beam that 
actually illuminates the reflector is determined. Then the representation of the 
beam on the dominium of the visibility diagram is identified. 

Replacing the point identified by the virtual source in the equation 
y = ax + b the equation of a straight line in the space of the parameters (a, b) is 
obtained, which together with the condition bo < b < bi represents a segment. 
Intersecting this segment with the areas that compose the visibility diagram of 
the reflector a series of sub-segments of the original segment is obtained, each 
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of which represents in the geometric space a new beam that will have incidence 
on a single reflector. Using this information a new level of the tree will be 
created. 

Figure 3 refers to the construction of a level of the tree, based on the 
reflections of the reflector u of Figure 2. Positioning the straight line S' in the 
space of the parameters (a, b), this intersects the reference area (-1 < b < 1) in 5 
sub-rays numbered from 1 to 5. The rays 1 and 3 are rays that are not reflected 
and proceed without limits, while the rays 2, 3 and 5 are reflected by the 
reflectors r2, ri and r3 respectively. The information relating to the position of 
the source that generates the reflection to that branch node, the geometric 
parameters that characterise the extension of the beam, and the information 
needed to identify the eventual reflector on which the beam will have incidence 
in continuing its path are associated to each branch node. In this manner the tree 
is created, and the process is repeated for each reflection of the rays, until the 
number of reflections reaches a preset number or the energy associate to the ray 
becomes negligible. 

The information present in the beam tree enable the generation of a filter 
(FIR) on the basis of which paths reach the location of the listener. This 
updating is rapid and can be carried out frequently. 

Through this process it is possible to obtain for each path its own direction 
of arrival (DOA), the length of the path and the list of the reflectors with which 
it has interacted along the path. This phase requires all the branches and the 
leaves of the tree to have been visited and has to be controlled if the receiver is 
present in the beam being considered. If the receiver is present, the direction of 
arrival (DOA) of the ray and the length of the path are calculated. Therefore, 
passing along the whole tree to the root the list of reflectors met is obtained. 

Since it is impossible to produce a flow for each DOA, angular intervals 
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have been grouped together to generate a limited number of audio flow, one per 
interval. Each of these angular intervals will be attributed to separate delay lines 
of the filter FIR, that are calculated taking into account all the paths and the 
relative attenuations. 

The filter calculated in this manner is capable, given an anechoic source 
(that is a source recorded in a room without reverberation), of obtaining the 
result of listening to that sound in a room represented by the geometry of the 
reflectors. 

The method allows immediate extension to the three-dimensional case. 
Also in the 3D case it is possible to carry out an affine transformation of the 
geometric space that enables the scansion interval to be preset, and to simplify 
the representation of the visibility function. As can be seen hereinafter, without 
losing generalities, a geometric transformation will be considered that brings the 
plane on which the reflector lies to coincide with the plane x=0 and brings the 
polygon to assume a preset canonic form. 

As already previously mentioned, the space of the parameters of the straight 
lines in three-dimensional environments is four-dimensional. Downstream from 
the geometric transformation of the polygonal reference reflector, the generic 
straight line of the space 3D can be expressed with the system of equations 
y = a x + b and z = c x + d, in which (a, b, c, d) identify the space of the 
parameters. In this space the representation of the potential visibility is still 
given by the set of the straight lines that have incidence on both the reflectors 
considered, which this time is a closed portion of the four dimensions space. For 
practicality, it is also possible to limit before the category of polygons that 
interest. A particularly advantageous choice is that of the triangles: this 
transformation does not remove generalities from the method, as it is always 
possible to represent a complex plane polygon by means of its triangulation; in 
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addition it is quite frequent, also in Computer Graphics, that the surfaces in the 
three-dimensional space are modelled with triangle nets. In the case of two 
triangular reflectors, the region of potential visibility can be expressed in the 
form of 6 disequation system, needed for the set of the straight lines having 
incidence on both the triangles. For example, expressing the generic reflector by 
means of the system of equations: 

x = a x t + b x s +c x 

y = a y t + b y s +Cy 

z = a z t + b z s +c 2 

t>0 

s>0 

t + s< 1 

The system that represents the generic straight line passing through the 
reference reflector, previously normalized by an affine transformation that 
makes its three vertexes coincide with the three points (0, 0, 0), (0, 1, 0) and 
(0, 0, 1), is represented by the system: 

y =ax+b 

z=cx+d 

b>0 

d>0 

b + d<l 

Replacing the first three equations of the first system, in the equations of 
the second system two equations are obtained, from which the parameters s and 
t can be derived in function of the variables a, b, c, d. Replacing the values of 
the parameters s and t in the disequations of the first system, and adding the 
disequations of the second system a system of 6 disequations is obtained in the 
parameters a, b, c, d, that delimit the area of potential visibility between the two 
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reflectors. 

In a similar manner to the two-dimensional case it is possible to decide the 
order of occlusion between overlapped regions of potential visibility, and thus it 
is possible to construct for each reflector a visibility diagram represented by a 
subdivision of the space of the parameters in regions of effective visibility, each 
of which is associated at the most to a reflector. 

The representation of the source and of the illuminated area is given by 
describing all the straight lines that intersect it with a limited beam, obtaining a 
system of two equations and some disequations in the four unknown quantities, 
that identifies a portion of a two-dimensional subspace. The subdivision of the 
beam will come about in this case intersecting the areas of the visibility diagram 
associated to the reflector that has generated the beam with the portion of plane 
that represents the beam reflected in the space (a, b, c, d). The method for the 
construction of the tree remains thus completely similar to the two-dimensional 
case. 

It can be noted that the choice of the transformation on the plane x = 0, both 
in the two-dimensional case and in the three-dimensional case, is not binding, 
but it has been introduced with the sole object of simplifying the exposition, 
placing the equations of the straight lines in an easier form. 

The method for determining the visibility between reflectors and the 
method for determining a beam tree, herein described, can be transcribed in a 
program, in a well-known manner by a technician of the sector, and can be 
memorised on any type of memory or support (floppy, CD) and/or can be made 
to be carried out by a computer. 
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